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ABSTRACT 

SGR 0418+5729 is a transient Soft Gamma-ray Repeater which underwent a major outburst in 
June 2009, during which the emission of short bursts was observed. Its properties appeared quite 
typical of other sources of the same class until long-term X-ray monitoring failed to detect any period 
derivative. The present upper limit on P implies that the surface dipole field is B p < 7.5 x 10 12 G 
(Rea ct al. 2010), well below those measured in other Soft Gamma-ray Repeaters (SGRs) and in the 
Anomalous X-ray Pulsars (AXPs), a group of similar sources. Both SGRs and AXPs are currently 
believed to be powered by ultra-magnetized neutron stars (magnetars, B p « 10 14 -10 15 G). SGR 
0418+5729 hardly seems to fit in such a picture. We show that the magneto-rotational properties 
of SGR 0418+5729 can be reproduced if this is an aged magnetar, m 1 Myr old, which experienced 
substantial field decay. The large initial toroidal component of the internal field required to match the 
observed properties of SGR 0418+5729 ensures that crustal fractures, and hence bursting activity, can 
still occur at present time. The thermal spectrum observed during the outburst decay is compatible 
with the predictions of a resonant compton scattering model (as in other SGRs/AXPs) if the field is 
low and the magnetospheric twist moderate. 
Subject headings: sources (individual): SGR 0418+5729 — stars: magnetic fields — stars: neutron 



1. INTRODUCTION 

Soft Gamma Repeaters (SGRs) and Anomalous X- 
ray Pulsars (AXPs) form a small, but rapidly expand- 
ing, class of isolated neutron star (NS) sources char- 
acterized by the emission of short (~ 0.1 s), energetic 
(« 10 40 ergs -1 ) bursts of X-rays. Their persistent lu- 
minosity (L x ~ 10 32 -10 36 ergs -1 in the ~ 0.5-10 keV 
range) is often variable, with flux enhancements up to 
several hundr eds during the outburs t phases of transient 
sources (e.g. IRea fc Esposito! 120111) . SGRs and AXPs 
share similar timing properties, with periods in a nar- 
row range (P ~ 2-12 s), and large period derivatives 
(P k, 10 -13 -10 -10 s s -1 ). In most of these sources 
the X-ray luminosity exceeds the rate of rotational en- 
ergy losses (E ?s 10 32 -10 35 ergs -1 ), and highly variable 
radio activity has been detected in three objects. X- 
ray spectra are often characterized by a thermal (BB; 
kT w 0.5 keV) plus a high-energy po wer-law (PL; pho - 
ton index T « 1.5-3) component (e .g. IRea et al.l 12 008). 
or by two thermal components (e.g. IGotthelf fc Halpernl 
120071: iTiengo. Esposito fe Mereghettil 120081 ). Both the 
spectral parameters and pul se profiles appear to vary 
with the source flux (see e.g. iWoods fc Thompsonll2006t 
IMereghettil 120081: IRea fc Esposito! 1201 11 for reviews on 
SGRs/AXPs properties) 6 . 
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The large values of the magnetic field derived from 
the standard dipole formula (5 x 10 13 G < B p < 
10 15 G), the lack of detected stellar companions and 
their large X-ray output as compared to E led to 
the suggestion that SGRs and AXPs are powered by 
a young (age « 10 3 -10 4 yr), ultra-magnetized neu- 
tron star, or magnetar (iDuncan. fc Thompson! Il992t 
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iThompson fc Duncan! [1993). Indeed, the magnetar sce- 
nario has been largely successful in explaining many of 
the observed p roperties of SGRs/AXPs, i ncluding those 
of the bursts (jThompson fc Duncan! 119951 ) and of the 
persistent emission (the twisted magnetosphere model , 
Thompson. Lvutikov fc Kulkarnill2002l: iZane et al.ll2009t 



Albano et al.ll201Ct and references therein). 



Despite SGRs and AXPs are far from being a homo- 
geneous class, in particular the inferred surface dipolar 
field spans nearly two orders of magnitude, their obser- 
vational behaviour is now commonly associated to that 
of (active) magnetars, to the point that often the terms 
SGR/AXP and magnetar are used as synonyms. This, 
actually, reflects the original definition of a magnetar 
as a NS which is powered b y its (large) magnetic field 
(Thompso n fc Duncan! I1993J ). In this respect, a super- 
strong magnetic field is not per se a sufficient condi- 
tion for triggering SGR/AXP-like activity, as testified by 
the existence of NS sources, for instance mo st of the so - 
called High-B radio pulsars (HBPSRs; e.g. lKaspil feOlO). 
and possibly s ome of the th ermally emitting isolated NSs 
(XDINSs; e.g. lTurollal[2009D . with surface magnetic fields 
comparable to those of SGRs/AXPs but having substan- 
tially different properties and not showing any burst- 
ing/outbursting activity over the ~ 10-20 yr time span 
during which they were observed. 

Over the last few years, the family of the mag- 
netar candidates grew with the addition of several 
new objects, most of them transient. X-ray burst- 
ing activity or peculiar radio emission similar to that 
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of SGRs/AXPs was detected in allege dly rotation- 
powe r ed pulsars, as PSR J18 46-0258 (jGavriil et al.l 
[20081: IKumar fc Safi-Harbl 120081) and PSR 1622-4950 
(jLevin et al.1 120101 k This suggests that magnetic en- 
ergy substantially contributes to their emission at certain 
stages. The magnetic field inferred from their rotation 
properties (B p ~ 0.5-3 x 10 14 G) is, in fact, close to 
that of SGRs/AXPs, contributing to the widespread be- 
lief that magnetar-like activity has to be associated to 
super-strong magnetic fields, typically higher than the 
quantum field B Q ~ 4.4 x 10 13 G. 

In this respect, t he recent discove ry of a low-field SGR, 
SGR 0418+5729 (jRea et al J 1201(1 . came as a surprise. 
SGR 0418+5729 was observed for the first time in June 
2009 when it entered an o utburst state during whic h 
X-ray bursts were detected (|van der Horst et aT1l2010f ). 
The enhanced flux level allowed for a measure of the 
source periodicity, P ~ 9.1 s, but, despite the source was 
then monitored in X-rays for ~ 500 d, no significant ev- 
idence for a period de rivative was found (jEsposito et al.l 
120101: iRea et al.ll2010L and references therein). The pub- 
lished upper limit is P < 6 x 10~ 15 ss" 1 , leading to an 
inferred dipole field B p < 7.5 x 10 12 G. 

The very low magnetic field of SGR 0418+5729 
(when compared to other SGRs/AXPs) raises a num- 
ber of questions as to if, and how, the observed 
phenomenology of this source can be accommodated 
within the magnetar picture. A crucial point is 
whether a NS with a surface field well below 10 13 G, 
can harbor an internal toroidal magnetic field strong 
enough to produce crustal displacements, which are 
believed to be responsible for the burstin g/outbu rsting 
episodes in SGRs/AXPs {T hompson & Duncan 1995; 
Thompson. Lvutikov fc Kulkarnil 120021 : iBeloborodov 



2009). In this paper we address this issue and dis- 



cuss both the spectral and timing properties of SGR 
0418+5729 in the framework of an aging magnetar. 

2. SGR 0418+5729 AS AN OLD MAGNETAR 

The suggestion that SGR 0418+572 9 may be an aged 
magnet ar was already put f orward by IRea et all (|2010l 
see also lEsposito et al] |2010) on the basis of the low per- 
sistent luminosity (most likely well below that observed 
during the outburst decay, Lx ~ 6x 10 31 (Z>/2kpc) 
ergs^ 1 ), weak bursting activity (only two faint bursts 
were recorded, van dcr Horst ct al. 2010), and large spin- 
down age (t c > 24 Mvr). IRea et al.l (J2010D estimated 
that an internal toroidal field B tor ~ 5 x 10 14 G is re- 
quired to power the persistent emission over a source 
lifetime ~ t c , and concluded that B tor can be still large 
enough to overcome the crustal yield. 

The main appeal of the "old magnetar" scenario is 
that it can offer an interpretation of the observed prop- 
erties of SGR 0418+5729 within an already established 
framework, validating the magnetar model also for (sur- 
face) field strengths quite far away from those of canoni- 
cal SGRs/AX Ps. Howev e r intr iguing, the considerations 
presented by IRea et ail (J2010D necessary rely on quite 
crude estimates. A more thorough investigation on the 
behaviour of evolved magnetars, in which a substantial 
decay of the magnetic field occurred, is definitely in order 
before claiming that SGR 0418+5729 is indeed powered 
by the last hiccups of an once ultra-magnetized NS. 



In the following we focus on three points which we 
deem central in order to test the "old magnetar" hy- 
pothesis. First, assuming that SGR 0418+5729 was 
born with a magnetar-like surface field, B p must have 
decayed by a factor k, 100 to match the current up- 
per limit. Roughly the same reduction is expected in 
the internal field. Although the latter can initially be 
w 10-100 times higher than B p (at least locally), one 
may wonder if at late times internal magnetic stresses 
are still strong enough to crack the crust. A second 
and related question is if realistic models of field de- 
cay in magnetars can account for the observed rota- 
tional properties (period and period derivative) of SGR 
0418+5729. This also directly bears to the true age of the 
source which is most probably much younger than the 
characteristic age, estimated assuming a non-decaying 
field. A final point concerns the persistent emission of 
SGR 0418+5729. XMM-Newton observations show evi- 
dence for a two component thermal spectrum (similar to 
those observed in other transient magnetars), although 
the presence of a non-therm al tail can not be excluded 
([Espositoet al.H2010L 1201 11 see 3!l4]) . In the magnetar 
model, spectra are expected to exhibit a power-law tail 
which originate s from resonant scattering in a twiste d 
magnetosphere ([Thompson. Lvutikov fc K ulkarni 2002). 
However, no calculations have been performed yet in the 
range of surface fields implied by SGR 0418+5729. 

2.1. Magneto-rotational Evolution 

A major issue in establishing the magnetic evolution of 
NSs (and of magnetars in particular) is that observations 
place very little, if any, constraint on the structure and 
strength of the internal magnetic field. While there 
are several indications that the large-scale, external 
field can be reasonably assumed to be dipolar, with a 
moderate amount of twist in magnetars, the different 
mechanisms proposed for the generation of the internal 
field in the earlier phases of the NS life (differential 
rotation, dynamo, magneto-rotational instability) most 
likel y give rise to both toroidal and poloi dal components 
(e.g. IGeppert. Kuker fe Page 2004, 2006, and references 
therein). The presence of a toroidal field, roughly in 
equipartition with the poloidal o ne, is also required by 
gener al stability arguments (e.g. iBraithwaite fc Sprufti 
120061 and references therein). A further complication 
comes from the present poor knowledge of where the 
internal field resides. The field can either permeate 
the entire star ("core" fields), or be mostly confined 
in the crust ("crustal" fields), depending on where its 
supporting (super)currents are located. The highly 
anisotropic surface temperature distribution required in 
some XDINSs has been taken as observational evidence 
in favor of a complex field geometry in the external 
layers (crust, envelope, atmosphere) of NSs, either in the 
form of strong crustal to roidal fields, multipolar p o loidal 
components, or both (fGeppert. Kuker fc Page! 120061: 



Perez- Azorfn. Miralles fc Pons! 120061: iZane fc Turollal 
20061) . 



The more general configuration for the internal field 
in a NS will then be that produced by the superposi- 
tion of cur rent systems in the cor e and the crust. As 
stressed by iPons fc Geppertl (|2007f) . the relative contri- 
bution of the core/crustal fields is likely different in dif- 
ferent types of NSs. In old radio pulsars, where no field 
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decay is observed, the long-lived core component may 
dominate, while a sizeable, more volatile crustal field 
is probably present in magnetars, for which substantial 
field decay over a timescale ~ 1 3 -10 5 yr is expected (e.g. 
Goldrcich fc Reisenegger 1992). 

A particularly important result 

(|Glarnpedakis. Jones fc Samuelssonl 120 111) is the lesser 
role that ambipolar diffusion plays in magnetar cores 
(after crystallization, the absence of convective motions 
already quenched ambipolar diffusion in the crust) 
on their active lifetimes, contradicting an assumption 
often made in the modelling of the flaring activity. 
Therefore, if the decay/evolution of the magnetic field 
is indeed the cause of magnetar activity, it is likely to 
take place outside the core and will be governed by 
Hall/Ohmic diffusion in the stellar crust. The relative 
importance of these two mechanism is strongly density- 
and temperature- dependent. Thus, any self-consistent 
study of the magnetic field evolution must be coupled 
to a detailed modelling of the neutron star thermal 
evolution, and conversely. Other mechanisms, e.g. flux 
expulsion from the superconducting core, due to the 
interaction between neutron vortices and magnetic flux 
tubes, are highly uncertain and very difficult to model. 
For these reasons, recent investigations of the magnetic 
field evolution in magnetars focused only on the crustal 
component of the field. 

The first attempts in t his d irection used a split ap- 
proach. iPons fc Geppertl ()2007f ) studied the evolution of 
the field by solving the complete induction equation in 
an isothermal crust, but assuming a prescribed time de- 
pendence for the temperature. They found that crustal 
magnetic fields in NSs suffer significant decay during the 
first « 10 6 yr and that the Hall drift, although inherently 
conservative (i.e. alone it can not dissipate magnetic en- 
ergy), plays an important role since it may reorganize 
the field from the larger to the smaller (spatial) scales 
where Ohmic dissipation proceeds faster. 

The cooling of magnetized NSs w ith field decay was 
investigated by Aguilcr a et al.1 ([2008) by adopting a sim- 
ple, analytical law for the time variation of the field 
which incorporates the main features of the Ohmic and 
Hall terms in the induction equation. The fully coupled 
ma gneto-thermal evolution of a NS w as finally addressed 
by IPons. Miralles fc Geppertl (|2009l) , including all real- 
istic microphysics. However, owing to numerical diffi- 
culties in treating the Hall term, their models include 
only Ohmic diffusion. This can be a limitation because, 
as they note, the Hall drift likely drastically affects the 
very early evolution of ultra-magnetized NSs with sur- 
face field B p > 10 15 G, and also that of "normal" NSs 
at late times (> 10 6 yr), when the temperature in the 
crust has dropped. On the other hand, for initial values 
of B p < 10 15 G, still well within the magnetar range, 
the effect of the Hall drift is expected to introduce at 
most quantitative changes (a somewhat faster dissipa- 
tion) with respect to the purely Ohmic picture. 

2.2. The Case of SGR 0418+5729 

To explore if, and to which extent, the magneto- 
thermal evolution of (initially) highly magnetic NSs can 
lead to objects with properties compatible with those of 
S GR 0418+5729, we perform ed som e runs using the code 
of IPons. Miralles fc Geppertl (|2009fl . We refer to Section 



2 in IPons. Miralles fc Geppertl (|2009l ) and Section 4 of 
lAguilera et al.l (|2008l ) for all details about the code and 
the microphysical input. We evol ved a lAMp* NS a ssum- 
ing the minimal cooling scenario ([Page et al.l l2004). with 
no exotic phases nor fast neutrino cooling processes, but 
including enhanced neutrino emission from the breaking 
and formation of neutron Cooper pairs in the NS core, as 
recent observations of the Cassiopeia A s upernova rem- 
nant seem to require (Pag e et al.l 120111 : iShternin et al.l 
1201 lh . The initial period was fixed at 10 ms and the 
initial dipole field to B p = 2.5 x 10 14 G. Note that the 
internal poloidal field is actually higher, with a maximum 
value B po i(t — 0) ~ 2.5 x 10 15 G in the inner crust. 

We considered three models with different values of 
the (maximum) internal toroidal field, B tor (t = 0) = 
0, 4 x 10 14 and 4 x 10 16 G, which turns out to be the 
crucial parameter, as shown in Figure [T] The four pan- 
els illustrate the evolution of luminosity, dipole field B p , 

period P and period derivative P. Indeed the prop- 
erties of SGR 0418+5729 are recovered in the case of 
B tor (t = 0) = 4 x 10 16 G and age ~ 1.5 x 10 6 yr. The 
main conclusions can be summarized as follows. 

i) The low quiescent luminosity is easily explained con- 
sidering that the object is relatively old: even a NS born 
as a bright, hot magnetar becomes cool and dim at this 
age. 

ii) On the other hand, the observed period constrains 
the dipolar field: the initial dipole field can not be much 
higher that the one considered here in order to prevent 
the star from spinning-down too fast and reach periods 
longer than that observed at present. Obviously there 
are other large uncertainties, such as the angle between 
rotation and magnetic axis, that may reduce the period 
(we assumed an orthogonal rotator here) . 

iii) Although the components of the initial internal field 
Btor(t = 0) can be varied to some extent, a quite large 
value is required. A large toroidal field, in fact, implies 
strong currents, which, in turn, produce more heating 
and higher temperatures. This drives a faster global field 
decay, which makes it possible to match the observed 
upper limit on P and B p . 

We stress that, while there are no stringent arguments 
against such large internal fields in the NS crust, their 
real occurrence in magnetars is an open issue. A possi- 
bility is that if the Hall drift becomes very important and 
it results in much faster dissipation, then one can obtain 
the same results starting with lower initial toroidal fields. 
Finally, we note that the Hall term is bound to become 
important again for objects like SGR 0418+5729 at late 
stages (> 1 Myr) as the star cools down and the con- 
ductivity increases by several orders of magnitude. No 
calculations are available in this regime but the possible 
occurrence of a second "Hall-active" phase could lead to 
enhanced bursting activity and rapid field decay. This 
may b e an indication that t he estimate of the bursting 
rate in iPerna fc Pons! (|2011l ) is a lower limit. 

2.3. Occurrence of Bursts 

Very recently IPerna fc Pons! (1201 lh use d the magnetic 
evolution code of IPons fc Geppertl (120071) together with 
the cooling models by IPons. Miralles fc Geppertl (2009) 
to compute the magnetic stress acting on the NS crust 
at different times. Their baseline model has B p (t = 0) = 
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that discussed in iPerna fc Pons! (|2011f) , and we estimate 
the typical lapse time between events for an object like 
SGR 0418+5729 is ~ 20-50 years. 



2.4. Persistent Emission 



Fig. 1. — From top left to bottom right, the evolution of the luminosity, surface dipole field, period and period derivative according 
to the model discussed in the text. The three cases refer to Btor(t = 0) = (solid lines), Btor(t = 0) = 4 X 10 14 G (dotted lines) and 
B tor (t = 0) = 4 X 10 16 G (dashed lines). 

8x 10 14 G and B tor (t = 0) = 10 15 G. They found that the 
occurrence of crustal fractures (and hence of bursts) is 
not restricted to the early NS life, during which the sur- 
face field is ultra-strong, but can extend to late phases 
(age w 10 5 -10 6 yr; see their figure 2). Both the energet- 
ics and the recurrence time of the events evolve as the 
star ages. For "old" magnetars about 50% crustal frac- 
tures release rs 10 41 erg and the waiting time between 
two successive events is ~ 1-10 yr. They also made a 
longer run with a model with B p (t = 0) = 2 x 10 14 G 
and B tor (t = 0) = 10 15 G, for which the event rate is 
about a factor 10 smaller. 

The model we considered in Sj2.ll as representative of 
SGR 0418+5729 has B p (t = 0) very close to this lat- 
ter configuration, while B tor (t = 0) is larger. The 
present (maximum) value of the internal toroidal field is 
~9x 10 14 G and, although we did not perform any de- 
tailed simulations, we argue that the bursting rate of our 
m odel, at its present ag e, is similar to the second model 
of IPerna fc Pon s (2011), because the internal configura- 
tion of the magnetic field is similar. It is important to 
notice that, despite the much larger initial toroidal fields 
of the model presented in this paper, this leads to faster 
decay and therefore similar values are reached when the 
NS is a million years old. Comparing both models at the 
estimated age of 1.5 Myr, the internal toroidal field of the 
model presented in this paper is only twice larger than 



In order to investigate the spectral properties of the 
persistent emission from SGR 0418+5729 and its time 
evolution, we analyzed eight Swift XRT and one XMM- 
Newton EPIC spectra (see Table SI of IRea et all [2010l 
for more details). Preliminary res ults for some o f these 
datas ets were already reported in lEsposito et al.l (|2010l 
120111) . Spectra for all the epochs were fitted simul- 
taneously using XSPEC v. 12.6, with the value of the 
column density Nh tied within the different observa- 
tions. Several one- and two-component models were 
tried, including a sing le blackbody, a one-dimensional 
(RCS.IReaet alj [2008h and a three-dimensional (NTZ, 
iZane et al.l I2009T) resonant scattering model, a double 
blackbody and a blackbody plus power-law. All single 
component models give rather poor fits. While for the 
RCS and NTZ models the residuals and the x 2 values 
were not acceptable, a single BB decomposition properly 
reproduces all the data but the XMM-Newton spectrum, 
which is the only responsible for the relatively large \ 2 of 

7 Each Swift dataset contains several individual observations in 
order to obtain good enough statistics. 
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the simultaneous fit. Since the highest-quality available 
spectrum argues against the source having (at least at 
the early stage of the outburst) a single thermal spec- 
trum, we decided to add a second component to the 
multi-instrument fit. 

A BB+BB and a BB+PL model with all parame- 
ters free (except 7V#, see above) provide acceptable fits 
of comparable quality (Xred = 1-15 for 601 dof and 
xl e d = 1-12 for 601 dof, respectively). However, we 
stress again that both these spectral representations con- 
tain a large number (32 in total) of free parameters, and 
those associated to the second component are not re- 
quired by the seven Swift observations but only by the 
XMM-Newton one. Moreover, despite on a statistical 
ground there is no reason to prefer the BB+BB over 
the BB+PL model, we note that in the latter : i) the 
spectral index T changes dramatically and in a totally 
erratic way from one observation to another, and ii) T 
can be as large as ~ 6, arguing against a power-law as a 
physically-motivated representation of the second spec- 
tral component. On the other hand, the values of the 
spectral parameters in the BB+BB model appear to be 
reasonable and their time evolution is monotonic (see be- 
low). 

Since in the BB+BB best fit model the temperature 
and normalization of the colder BB component appear 
not to vary sensibly in time (again, possibly because 
they are poorly constrained by the Swift observations), 
we performed a fit with these two parameters tied across 
the various datasets. This resulted in a similarly good 
fit (Xred = 1-18 for 617 dof) and has the advantage to 
contain 16 degrees of freedom less. In the case of the 
BB+PL model the goodness-of-the-fit worsens consider- 
ably {xfed = 1-42 for 617 dof) by requiring that the pa- 
rameters of the (single) BB are the same at the different 
epochs. 

For these reasons, in the following we take the BB+BB 
model (with the colder BB constant in time) as the most 
likely representation of the data, and discuss the ensuing 
implications in framework of an evolved magnetar. It is 
worth mentioning that the failure of the resonant scat- 
tering models to fit the data may be due to the fact that 
both RCS and NTZ were originally developed for much 
higher fields than that likely present in SGR 0418+5729 
(the NTZ version used here assumes B p — 10 14 G). A 
more detailed spectral analysis will be the subject of a 
forthcoming paper (Rea et al. in preparation). 

The picture which emerges from the spectral analy- 
sis is that of thermal emission from two regions on the 
star surface, a cold one, with more or less constant size 
and temperature (R c ~ 0.75 km for a fiducial distance 
D = 2 kpc and T c ~ 0.31 keV), and a hot one, which 
shrinks during the outburst decay. The evolution of the 
temperature and size of the two components is shown 
in figure [2] The temperature of the hot region is more 
or less constant at fcl\ ~ 0.93 and its area changes from 
~ 0.2 to ~ 0.03 times that of the cold region (Rh ~ 0.15- 
0.30 km, again for D = 2 kpc). The overall behavior is 
quite reminiscent of those seen in other transient mag- 
netar sources, notably the A XPs XTE J1810-197 and 
CXOU J164710. 2-455216 (e.g. lAlbano et all \WHJ and 
references therein). 

Within the magnetar model this is interpreted as due 
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Fig. 2. — The time evolution of the temperature and emitting ra- 
dius of the two BB components in the spectrum of SGR, 0418+5729; 
a source distance of 2 kpc is assumed. Diamonds refer to the hot 
and filled circles to the cold component. The solid line shows the 
ratio of the emitting areas, A^^/A^i^. Time is counted in days 
from the outburst onset, on 2009-06-05 20:38:24.000 UTC (MJD 
54987.862). 

to the sudden development of a twist in the external 
magnetic field which then progressively decays. The 
twist likely affects only a limited bundle of (closed) field 
lines and the charges flowing along the current-carrying 
bundle heat the surface layers as they impact upon the 
star. When the magneto sphere untwists, th e size of the 
heated region decreases (jBeloborodovl I2009J) . This pic- 
ture is compatible with the results we obtained for SGR 
0418+5729 with the BB+BB model assuming that the 
heated region corresponds to the area emitting the hotter 
BB. This is superimposed (or close) to a cooler, larger cap 
which is responsible for the emission of the softer BB. It 
is interesting to note that the analysis of the pulse profiles 
during the first stages of the outburst supports this view. 
The double-peaked pulse profile of SGR 0418+5729 sug- 
gests, in fact, that the surface thermal map of the star 
comprises two warm caps, only one of which was involved 
in the heating process (Esposito et a l. 2010). The pre- 
dicted characteristic time for the outburst evolution is 
« 5($/lO 9 V)^ 1 (B/lO 14 G)A0(A/lO 12 cm 2 ) yr, where $ 
is the discharge voltage and A is the area of th e surface 
region involved by the twist (|Beroborodovll2009l) . Taking 
B~5x 10 12 G, A0 - 0.4 rad and A ~ 10 11 cm 2 , we 
get for the characteristic time » 0.02($/10 9 V) _1 y r - A 
low discharge voltage, $ « 10 8 V, is then required to 
obtain a decay tim e ~ a few mon t hs. W e warn that, 
as already noted bv lEsposito et all (|2010f ). the luminos- 
ity produced by ohmic dissipation appears to be too low 
to reproduce that observed at the beginning of the out- 
burst, ~ 10 34 (D/2kpc) 2 ergs -1 . However, if the twist 
affected a region different from a polar cap (e.g. a ring 
confined between two values Q\ and 6*2 of the magnetic 
colatitude) the value of the luminosity can be higher 8 . 
Alternatively, other heating mechanisms may be operat- 

8 We thank A. Beloborodov for bringing this point to our atten- 
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ing, e. g. the release of magnetic energy in the star outer 
layers ( Lvubarskv. Eichler fc Thomp son 20Q2J). 

As discussed in Section [2731 the internal field of SGR 
0418+5729 can be still large enough to produce crustal 
displacements so that magnetic helicity is transferred 
to the external field, twisting up the magnetosphere. 
The appearance of a twist is usually accompanied by 
the formation of a high-energy spectral tail, due to res- 
onant cyclotron up-scattering of thermal surface pho- 
tons, which is, however, not unambiguously detected in 
SGR 0418+5729. In order to investigate the properties 
of resonant cyclotron scattering spectra in the low-field 
regime (B p < 5 x 10 13 G), we run a series of 3D Monte- 
carlo simulations, using the rel ativistic transport code of 
iNobili. Turolla fc Zand (|2008allbT ). to which we refer for 
all details. We considered four values of the (polar) sur- 
face field {B p = 10 12 , 5 x 10 12 , 10 13 , 5 x 10 13 G), two val- 
ues of the seed photon temperature (kT = 0.3, 0.9 keV) 
and several values of the twist angle 9 in the range 
0.1 rad < A</> < 1.2 rad. The electron temperature and 
bulk velocity were fixed to kT e i = 10 keV and v/c = 0.5 
in all cases. A different choice of these parameters pro- 
duce similar results provided that the scattering particles 
are mildly relativistic, as indeed required to reproduce 
the observed 1-10 keV spectra o f SGRs/AXPs (se e e.g. 
INobili. Turolla fc Zand l2008allbl ; iZane et all [200l and 
references therein) . Results are summarized in figure [3] 
which shows the photon index of the non-thermal tail 
(computed in the 6-8 keV range) as a function of the 
twist angle for the different values of B p . The (average) 
index (in the same energy range) of the blackbody spec- 
trum is marked by a dashed horizontal line: when the 
photon index approaches the line the spectrum becomes 
indistinguishable from a blackbody and no tail is present. 
As it can be seen, while for kT = 0.3 keV a non-thermal 
tail below 10 keV appears for all the values of the twist, 
unless B p = 10 12 G and A<j) < 0.3 rad, the up-scattering 
of seed photons associated to the hotter component only 
produces a tail if Acj) > 0.5 rad and B p > 10 12 G. We 
stress that here we are considering photon energies below 
10 keV, so the lack of a non-thermal spectral component 
for kT — 0.9 keV only reflects the fact that now resonant 
comptonization tends to move photons at energies higher 
than 10 keV. A tail, in fact, may be present above 10 keV 
also if it does not show up below 10 keV. 

Although the resonant scattering spectrum produced 
by the reprocessing of soft photons coming from two NS 
surface regions at different temperature is not exactly 
given by the superpo sition of the two ind ividual models 
(see the discussion in lAlbano et al. 2010), wc adopt this 
approach to get some insight into the spectral proper- 
ties of SGR 0418+5729. If the observed spectrum of the 
source is best modelled in terms of the superposition of 
two blackbodies with kT ~ 0.3, 0.9 keV the twist an- 
gle must be < 0.5 rad for B p > 10 12 G not to produce 
a power-law tail in the hot component (see figure |3]). 
This, however, is only a necessary condition because a 
PL tail may still appear in the cold component. The 
emergence of such a power-law is related to the rela- 
tive magnitude of the hot and cold components. The 
total spectrum resulting from the superposition of the 

9 Here the magnetosphere is assumed to be globally twisted. 



models with kT ~ 0.3 and 0.9 keV is shown in figure 
|4] for Bp — 5x 10 12 G and two values of the ratio of 
the emitting areas, A co id/Ah ot — 15, 30, typical of those 
measured during the evolution of SGR 0418+5729. We 
remark that the largest area ratio (= 30) corresponds to 
the most unfavorable case: if the tail does not appear 
now it is not present for smaller values of the area ra- 
tio, when the cold component contributes less. As it is 
seen, the total spectrum is very close to the superposition 
of two blackbodies, with no high-energy tail. The same 
result holds for different values of the magnetic field, pro- 
vided that A</> < 0.5 rad, and for even larger values of 
the twist if the field is as low as 10 12 G. We conclude 
that the strong evidence that SGR 0418+5729 exhibits a 
X-ray spectrum dominated by surface thermal emission 
below 10 keV is not in contradiction with the predictions 
of RCS models. 

3. DISCUSSION 

SGR 0418+5729 is the first "low-field" soft gamma re- 
peater/anomalous X-ray pulsar ever discovered. Even 
if other sources of the same class with a relatively 
weak i?-field were al ready known (e.g. IE 2259+586, 
iGavriil fc Kaspill2002D . the case of SGR 0418+5729 is ex- 
treme and stimulated considerable interest. With a sur- 
face dipole field < 7.5 x 10 12 G, SGR 0418+5729 seems 
to challenge an interpretation in terms of the "conven- 
tional" magnetar model, and not just because its mag- 
netic field is below the quantum critical field, a quite 
irrelevant fact per se. 

The present upper limit on the surface field in SGR 
0418+5729 is based on the spin-down measure, which 
traces the dipole component. It is actually possible that 
the external magnetic field in NSs, and in SGRs/AXPs 
in particular, is more complex than a simple dipole. 
Higher order multipoles can substantially contribute to 
the field near the surface and, because they fall off more 
rapidly with radius, induce negligible spin-down. SGR 
0418+5729 may, then, possess a much higher surface B- 
ficld than indicated by its spin-down rate. 

A high surface field in the form of multipolar com- 
ponents hints toward the presence a large internal field 
which can produce crustal motions and bursting activity, 
making SGR 0418+5729 not dissimilar (a part from the 
magnetic field topology) from other magnetar sources. 
It does not, however, explain the rotational properties of 
the source. With a dipole field below 10 13 G, in fact, 
it would be impossible to slow down the star to the ob- 
served 9.1 s period in less than ~ 24 Myr, a time much 
longer than the estimated age of other SGRs/AXPs, 
unless SGR 0418+5729 had an exceptionally long pe- 
riod at birth. A possible sol ution was suggested by 
lAlpar. Ertan fc Caliskanl (J2011D who have shown that a 
if SGR 0418+5729 was born with a period > 70 ms and 
a low dipolar field (B w 10 12 G), the torque exerted 
by a fallback disk can spin down the star to the present 
period in > 10 5 yr, if P has to be below the observed up- 
per limit 10 . The survival of multipolar field components 
sw 100 times stronger than the dipole over such a time 

10 A somehow similar scenario in which SGRs/AXPs are NSs 
with a low dipole field and super-strong multipolar components 
po wered by accret i on f rom a fallback disk was recently proposed 
bv llrumper et al.l (IMToT ). 



Is SGR 0418+5729 indeed a waning magnetar 
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Fig. 3. — Left panel. The photon index T vs. the twist angle A(f> for B p = 10 G (filled circles), 5 X 10 G (open circles), 10 3 G 
(squares) and 5 X 10 13 G (diamonds); the seed photon temperature is kT = 0.3 keV. The dashed horizontal line marks the value of the 
(average) index of the blackbody spectrum at the same temperature and in the same energy range (see text). Right panel. Same for 
kT = 0.9 keV. 



(Lx ~ 6 x 10 31 ergs 1 for a distance of 2 kpc; Re a et all 
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Fig. 4. — The spectra obtained from the superposition of two 
resonant scattering models with kT = 0.3, 0.9 keV, A</> = 0.4 rad, 
B p = 5 X 10 12 G. The cold component is shown in red, the hot in 
blue and the sum in green. Solid (dashed) lines are for an emitting 
area ratio A co i d /A hot = 15 (30). 



span may, however, be an issue in the light of the known 
evolution of the dipolar field in magnetars (see ij2.1|) . 

In this paper we explored a different possibility, i.e. 
that SGR 0418+5729 is an old neutron star born with 
a super-strong magnetic field which experienced field 
decay over a time « 10 6 yr. Our results show that 
magneto-dipolar braking can effectively spin down the 
star to a period ~ 10 s with P < 10~ 15 ss" 1 provided 
that the initial internal toroidal field is large enough, 
Btor(t = 0) > 10 16 G. At the same time the initial ex- 
ternal dipole field has to be < 2-3 x 10 14 G in order 
to prevent the NS from spinning down too fast. In par- 
ticular, the model with B tor (t = 0) = 4 x 10 16 G and 
B p (t = 0) = 2.5 x 10 14 G gives B p ~2x 10 12 G, P ~ 9 s 
and P~4x 10~ 16 ss" 1 at an age of ~ 1.5 Myr, in very 
good agreement with the current observational picture of 
SGR 0418+5729. The predicted quiescent luminosity at 
the same age is Lx ~ 10 31 ergs -1 , again in agreement 
with the current luminosity of the source, the faintest 
measured so far and possibly close to the quiescent value 



I2010f ). We note in this respect that the luminosity drops 
very quickly for ages > 10 6 yr (top left panel of figure [TJ, 
so the present value could well be below ~ 10 31 ergs -1 if 
the age is only slightly above 1.5 Myr. A somehow longer 
age would have negligible impact on the predicted P, P 
and B p which already reached a nearly constant value 
(see again figure []}. We also note that the measured 
flux in the ~ 0.5-10 keV band may not be representa- 
tive of the bolometric luminosity for an old object like 
SGR 0418+5729. If the NS surface temperature is < 10 6 
K, in fact, the quiescent thermal emission would be too 
soft (and absorbed) to be clearly detectable. It is then 
possible that the observed X-ray flux only constrains the 
"outburst" luminosity and the genuine quiescent emis- 
sion may go undetected even if it is w 10 31 ergs" 1 . 

Given the large internal toroidal field at birth, our fidu- 
cial model for SGR 0418+5729 retains the capability to 
induce crustal fractures, and hence to produce bursts, 
even at quite late times. A comparison of our model with 
one of the cases investigated in detail by iPerna fc Pons! 
(J2011I ) at a comparable age (~ 1.5 Myr), indicates that 
the recurrence time between bursting/outbursting events 
for an object like SGR 0418+5729 is a few tens of years. 
Crustal displacements are accompanied by a twisting up 
of the magnetosphere and the formation of a high-energy 
spectral tail, which is observed in most SGRs/AXPs. We 
have shown, however, that a source with a spectral dis- 
tribution consistent with a double blackbody, as follows 
from the analysis of a series of X-ray observations taken 
during the outburst decay of SGR 0418+5729, can be 
modelled as well by a resonant cyclotron scattering model 
if the twist is moderate (twist angle < 0.5 rad) and the 
surface field is low, < 5x 10 12 G, as predicted by our 
evolutionary calculations. 

It has been already noticed that the SGRs/AXPs 
which exhibit the larger flux variations seem to b e those 
with the lower dipole fields (Esposito et al. 2009). SGR 
0418+5729 provides a further case for such a correlation. 
In the magnetar scenario the occurrence and energetics 
of outbursts are dictated by the internal field and the 
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peak luminosity L max depends on the (maximum) en- 
ergy stored locally prior to the event. L max is not much 
sensitive to the external dipole field. If the energy re- 
leased in an event is roughly similar in all sources (as 
in the case the mechanism is gated e.g. by the crustal 
yield), the ratio between the peak and persistent fluxes 
has to be much higher in the less active, old objects than 
in active, young ones. This is because L max is similar for 
all of them, but the quiescent luminosity is much lower 
for old sources. In the latter, field decay had the time to 
reduce the dipole field to rather low values, so old, low- 
field objects (but with a sufficiently high internal field) 
appear as transient sources. It is intriguing that SGR 
0418+5729 is the source with the lowest field known and 



at the same time a most extreme transient, with a peak- 
to-persistcnt flux ratio ~ 1000. 
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